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We present here the microwave surface impedance of a high purity crystal of YBa^CuzO^m mea- 
sured at 5 frequencies between 1 and 75 GHz. This data set reveals the main features of the 
conductivity spectrum of the thermally excited quasiparticles in the superconducting state. Below 
20 K there is a regime of extremely long quasiparticle lifetimes, due to both the collapse of inelastic 
scattering below T c and the very weak impurity scattering in the high purity BaZrO^-grown crystal 
used in this study. Above 20 K, the scattering increases dramatically, initially at least as fast as T 4 . 
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I. INTRODUCTION 

Over the past few years measurements of electrody- 
namic properties at microwave frequencies have proven 
to be a particularly fruitful technique for studying the 
superconducting state of the high temperature supercon- 
ductors. A key strength of the technique is that mea- 
surements of the real and imaginary part of the surface 
impedance Z s (lu,T) provide complementary information 
on two aspects of the superconducting state; the super- 
fluid density and the low energy excitations out of the 
condensate. 

Measurements of the imaginary part of the surface 
impedance X S (T) provide a direct measurement of the 
penetration depth A(T) and thus give a direct measure 
of the temperature dependence of the superfluid density. 
Such measurements have produced a wealth of informa- 
tion on the superconducting state. The widespread ob- 
servation of a linear temperature dependence of A(T) at 
low T in many of the superconducting cuprates (l]-|| has 
been a key piece of evidence suggesting nodes in the en- 
ergy gap in these materials. Near T c the temperature de- 
pendence of the superfluid density has provided evidence 
of 3DXY critical fluctuations over a wide temperature 
range in Y Ba-iCuzO-is ||. 

Measurements of the real part of the surface impedance 
R S (T), when combined with the measurements of A(T), 
can be used to determine the real part of the microwave 
conductivity cti(T), which is essentially electromagnetic 
absorbtion by quasiparticles excited out of the conden- 
sate (either thermally excited quasiparticles or excita- 
tions created by the absorbtion of photons). Early mea- 
surements of R S {T) at a few GHz exhibited a broad peak 
below T Cl caused by a very large peak in a\ (T) @ at these 
low frequencies. This peak was also observed at THz fre- 
quencies by Nuss et al. |^| , and was attributed to a com- 
petition between two temperature dependences; the over- 
all decrease with temperature of the number of thermally 
excited quasiparticles competing with a rapid increase 
below T c of the transport lifetime of these quasiparti- 
cles. This rapid increase in lifetime has been interpreted 
as a collapse of the inelastic scattering processes respon- 
sible for the large normal state resistivity of the high 



temperature superconductors and is an effect that is ob- 
served in thermal conductivity measurements |)| as well 
as in other electromagnetic absorption measurements at 
microwave [pj,p~0 16 , far infrared JlTj] and THz frequen- 
cies ^ P^P^i - The rapid increase in quasiparticle life- 
time is well established by these measurements and the 
very long quasiparticle mean free paths resulting from 
this have been corroborated further by the thermal Hall 
effect measurements of Krishana et al. |20 . However, 
obtaining a quantitative determination of the scattering 
lifetime and the details of its temperature dependence 
has been hampered by the need to use models to inter- 
pret the existing microwave data. The problem has been 
that although the step from the R S {T) and A(T) to the 
conductivity a\ (T) is a matter of straightforward super- 
conductor electrodynamics when in the local limit, the 
extraction of a scattering time from ai(T) suffered from 
a dependence on an assumed model for the shape of the 
quasiparticle conductivity spectrum ai(u>). 

In this paper we present measurements at 5 different 
microwave frequencies, giving enough spectroscopic de- 
tail over a wide enough frequency range to produce a 
rather complete picture of the evolution of ai(u>,T) in 
the superconducting state. These results support the 
early model based on an ansatz that the thermally ex- 
cited quasiparticles have a nearly Drude-shaped conduc- 
tivity spectrum and now provide a much clearer measure- 
ment of the temperature dependence of the quasiparticle 
scattering rate below T c . We find that this scattering 
rate becomes extremely small and appears to become es- 
sentially temperature independent below about 20 K. At 
higher temperatures the scattering rate increases rapidly, 
initially at least as fast as T 4 . In the following section 
we will describe the techniques used to produce this in- 
formation and in Section EH we will present the results 
along with the details of the extraction of the microwave 
conductivity from the surface impedance measurements. 
In Section IV the microwave conductivity spectra will 
be presented along with the results of fits that generate 
the temperature dependence of the quasiparticle lifetime. 
In Section ^ we will point out some of the implications 
of these results and in particular compare them to the 
present literature on conductivity in a ci-wave supercon- 



1 



ductor. 



II. EXPERIMENTAL TECHNIQUES 

The property being directly probed in a microwave 
measurement on a superconductor is the complex sur- 
face impedance Z S (T) = R S {T) +iX s (T). One simpli- 
fying feature of the high temperature superconductors 
is that the superfluid response falls in the limit of local 
electrodynamics, so that R S (T) and X S (T) are related to 
the complex conductivity a (to , T) — o\ (to , T) — 102 (to , T) 
in a straightforward way. This simple limit arises from 
the very small coherence length in these materials, which 
guarantees that ( « A. A particularly simple regime of 
the local limit occurs when a% 3> o\ below T c and at low 
frequency, in which case one obtains the relations 



R s (T) = ^cj 2 X 3 (T)a 1 (cj,T) 
X S (T) = /xowA(T) 



(2.1) 



A more detailed discussion of the electrodynamics and 
the extraction of <ti(uj,T) from Z S (T) will be presented 
in Section III , but the equations above are useful for quick 
estimates and for understanding the main features of the 
microwave properties of these superconductors. 

The expression for R S (T) in Eq. 2.1 embodies what 
is most difficult about the microwave measurements. It 
contains both o\ and a term of the form w 2 A 3 , and physi- 
cally can be interpreted as the microwave absorption pro- 
cesses (cti) occuring within the rather shallow depth into 
which the microwaves penetrate (hence the term w 2 A 3 ). 
This screening length set by the superfluid makes R s ex- 
tremely small below T c , such that for small single crys- 
tals, one is forced to employ cavity perturbation in very 
high Q microwave resonators. This is most often achieved 
by the use of microwave cavities made from conventional 
superconductors cooled to low temperature. This fixed 
frequency type of measurement must then be performed 
with several different resonators if one is to build up a 
reasonably complete picture of the microwave conductiv- 
ity spectrum ai(u)). 

The data presented here involve measurements with 5 
different resonators spanning a range from 1 to 75 GHz 
and utilizing a number of variations on the basic method 
of cavity perturbation. The common feature is that all of 
the measurements have been performed on the same sam- 
ple, a thin plate of Y Ba2Cu^O(,.g^ oriented such that the 
microwave magnetic field of each cavity lay in the plane 
of the plate (H r f || b). This geometry has the advantage 
that demagnetization factors are quite small, so the sur- 
face current distributions are fairly uniform and are sim- 
ilar in all of the measurements, making comparison from 
frequency to frequency quite reliable. The disadvantage 
of this geometry is that although it mainly measures the 
surface impedance for currents running across the crys- 
tal face in the a direction, there is some admixture of the 



c-axis surface impedance coming from currents running 
down the thin edge of the crystal. However, we have 
previously shown that these effects are small for a thin 
crystal, because the temperature dependences of R S (T) 
and X. (T ) arc quite weak in the c-direction, except near 
T c |^l|,^2| • The absolute value of R s (T) for the c-direction 
is also relatively low J2^| . 

The measurements at 1 GHz were performed in a loop- 
ap resonator initially designed for measurement of A(T) 
y. Like most of the resonators used in this study, the 
loop-gap is plated with a Pb:Sn alloy which is super- 
conducting below 7 Kelvin and has very low microwave 
loss at 1.2 Kelvin. In the case of the 1GHz loop-gap, 
the cavity Q can be as high as 4 x 10 6 at low tem- 
perature. Another common feature of all of the mea- 
surements is that the sample is mounted on a thin sap- 
phire plate with a tiny amount of silicone grease, with 
the thermometry and sample heater located at the other 
end of the sapphire plate, outside of, and thermally iso- 
lated from the resonator. In this way the resonator can 
be held fixed at the regulated 4 ife bath temperature 
while the sample temperature is varied. A unique fea- 
ture of the 1 GHz loop-gap system is that the sample is 
held fixed in the resonator, with the sapphire plate and 
thermometry stage supported by a thin walled quartz 
tube which sustains the temperature gradient between 
the i He bath and the sample. This means that the sam- 
ple cannot be removed from the resonator during the 
measurements, but we find this restriction is necessary 
for the high precision measurements of X(T) which rely 
on the sample being held rigidly in a position that does 
not vary when the temperature is changed. Without 
this type of construction, motion of the sample and sap- 
phire plate in the fields of the resonator can give rise 
to experimental artifacts. This is because A(T) is de- 
termined from very small changes in the resonator fre- 
quency as the sample temperature is changed. This tech- 
nique gives highly precise and reliable temperature de- 
pendences relative to the base temperature, but is limited 
to measuring differences; AA(T) = A(T) - A(1.2 K) and 
AR S (T) = R S (T) - R s (1.2 K). In this paper we will not 
be focussing on the low temperature limit R S (T — > 0), 
which is a very small microwave loss in comparison to 
the loss associated with thermally excited quasiparticles. 
On the other hand, we do need to know X(1.2K) in or- 
der to extract microwave conductivities, so we obtain this 
value from infrared measurements p3fl . 

The 2 GHz measurements have been performed in a Nb 
split-ring resonator that will be described in more detail 
elsewhere 24 1 . Its geometry is related to the 1 GHz loop- 
gap resonator, with the key operational difference being 
that the sample can be pulled in and out of the resonator 
to measure the Q and resonant frequency of the unloaded 
cavity. Measurements at 13.3, 22.7, and 75.1 GHz have 
been performed in the axial microwave magnetic fields 
of the TEqu modes of three right-circular cylindrical 
cavities. In these higher frequency cavity perturbation 
measurements, the sample can also be withdrawn from 
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the resonator. This allows one to determine R S {T) from 
R S (T) <x 1/Qs — 1/Qo where Qo is the Q of the empty 
cavity and Q s is the Q with the sample inserted. A small 
correction for other sources of loss is made by measur- 
ing the sapphire holder and grease without the sample; 
the calibration of the absolute value of the surface resis- 
tance is obtained by measuring a Pb:Sn reference sample 
whose normal state surface resistance is governed by the 
classical skin effect. 

The sample used for all of these measurements was 
a single crystal of Y Ba2Cu^OQ^ grown by a flux 
growth technique using BaZrO^ crucibles ^5|. The 
BaZrO^, which does not corrode during the crystal 
growth, yields crystals of much higher purity [^5|,|6| (bet- 
ter than 99.99%) than those grown in more commonly 
used crucibles such as yttria-stabilized zirconia and alu- 
mina. The sample used for this was detwinned at 250 
C under uniaxial stress and was subsequently annealed 
at (350 C) for 50 days to produce a sample with nearly 
filled chain oxygen sites. This produces a sample with a 
slightly lower T c (88.7 K) than the maximum obtained 
near an oxygen concentration of 6.91 p5[ , but provides 
particularly defect-free samples without the oxygen va- 
cancy clustering discovered by Erb et al. |27| . The dimen- 
sions of the sample were initially 2x1x0.02 mm 3 for the 1 
GHz measurements, which have been reported elsewhere 
p8j. Subsequent measurements in the higher frequency 
resonators were performed on smaller pieces cleaved from 
the original crystal. 

III. EXPERIMENTAL RESULTS AND ANALYSIS 

Fig. [l| shows the microwave surface resistance of 
Y Ba2C'u30e.gg3 for currents running in the a direction. 
At all 5 of the frequencies shown in the figure, the rapid 
drop in R S (T) below T c is due to the onset of screen- 
ing by the superfluid, with the overall magnitude of the 
drop depending strongly on frequency as expected from 
the term in R S (T) that varies as ui 2 A 3 (T). For easier 
comparison of the different frequencies, it is convenient 
to plot the low loss data in the superconducting state as 
Rs/oj 2 versus temperature, as shown in Fig. ^. One fea- 
ture of these figures is that they demonstrate that above 
65 K the loss scales quite closely as lu 2 , indicating that 
<7i is frequency independent above 65 K in the microwave 
range. Alternatively, if one expects u> 2 scaling, then this 
indicates that the comparison between the frequencies is 
working extremely well, even though the data is gener- 
ated in 5 different resonators with separate calibrations of 
the absolute value in each experiment. The curves in Fig. 
^ show a number of qualitative features that have already 
been observed in earlier microwave measurements of this 
material. R S (T) exhibits a broad peak which shifts up in 
temperature and diminishes in size as the measurement 
frequency is increased. Also noteworthy is that we con- 
tinue to observe nearly linear behaviour of R s (T) at low 
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FIG. 1. The surface resistance R S (T) for currents running 
in the a direction of a high purity crystal of YBa2Cu3.Ofi.gg3 
The measurements range over 5 different frequencies from 1 
to 75 GHz. 

temperature and low frequency. 

The broad peak in R S (T) was originally attributed to 
a quasiparticle scattering time r that increases rapidly 
with decreasing temperature below T c and competes with 
a density of these thermally excited qusiparticles that 
decreases with temperature 0. At intermediate tem- 
peratures from T c down to about 30 K the quasiparticle 
lifetime increases much more rapidly than the density de- 
creases, causing <7i(T) (and thus R S (T)) to rise with de- 
creasing temperature. For the early measurements near 
2 GHz it was suggested that the quasiparticle scatter- 
ing time reaches a limiting value near 30 K, possibly due 
to impurity scattering, at which point the quasiparticle 
density takes over and causes R s (T) to fall again with 
further decreases in temperature fl29| . The speculation 
that impurities are involved in the turnover at 30 K was 
partly checked by studying samples doped with Ni and Zn 
. These doping studies showed either a smaller peak 
shifted to higher temperature or no peak at all, consis- 
tent with the quasiparticle scattering rate running into 
an impurity limit at higher temperature, even when only 
0.15% Zn was added to the sample. This sensitivity to 
such low levels of impurities raises a serious concern over 
these earlier crystals grown in yttria-stabilized zirconia 
crucibles, because during crystal growth the residual im- 
purity level due to uptake of material from the corroding 
crucible reaches the 0.1% level. The results shown here 
on a new higher purity sample confirm the original spec- 
ulation that it is this residual impurity scattering that 
limits the increase in quasiparticle lifetime, even in quite 
pure crystals. The low frequency surface resistance of the 
new crystal rises considerably higher above the minimum 
near 70 K than was observed in earlier measurements 
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FIG. 2. The same measurements of surface resistance as 
shown in Fig. 1, but with the behaviour below T c emphasized 
by dividing out a frequency-squared dependence associated 
with superfluid screening. 

of samples at 2 and 4 GHZ and reaches its maximum 
at a lower temperature p9fl . Both of these quantitative 
changes are consistent with the new samples simply hav- 
ing lower impurity scattering and so reaching a higher 
quasiparticlc lifetime limit at a lower temperature than 
was seen in the yttria-stabilized zirconia grown crystals. 

The frequency dependence of the peak in R s (T) pro- 
vides further evidence of the rapid increase in the quasi- 
particle lifetime. In much higher frequency measure- 
ments on thin films, Nuss et al. || were the first to ob- 
serve a similar broad peak in the THz range that shifted 
up in temperature and decreased in size at higher fre- 
quencies. They pointed out that this could be accounted 
for by relaxation effects. If the quasiparticle lifetime in- 
creases sufficiently that lot > 1, then the sample enters 
a relaxation regime where <J\{T) falls as r continues to 
increase. So for high frequency measurements the tem- 
perature at which R S (T) reaches its peak roughly indi- 
cates where lot — 1 at each measurement frequency. The 
fact that we already observe these relaxation effects at 
13 GHz indicates that the scattering time r grows to ex- 
tremely large values in the new high purity samples. 

To better understand the data it is desireable to ex- 
tract cfi{lo,T) from these measurements of R s (lo,T). At 
the lowest frequency of 1 GHz, where lot <C 1 and we 
have simultaneous measurements of R S (T) and A(T), it 
is straightforward to extract o\ (T) with only an assump- 
tion that the electrodynamics are local. The 1 GHz mea- 
surements of A(T) used in this analysis are shown in Fig. 
|, plotted as 1/A 2 (T). The screening by the superfluid 
follows the local London model if a superconductor is in 
the limit (<Cl Strictly speaking, this local limit is a bit 
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FIG. 3. The temperature dependence of the o axis pene- 
tration depth of Y Ba2Cu2.Oe.95 plotted as 1/A 2 (T), which is 
a measure of the superfluid density via n B e 2 /m* = l/(/^oA 2 ). 

more complicated for the case of a superconductor with 
nodes in the energy gap because the coherence length 
is then /c-dependent and becomes large in the node di- 
rections. However, the consequences of this type of non- 
locality for the electrodynamics would only be observable 
at low temperatures and the effect would be small in the 
measurement geometry used here pj] . A much more se- 
rious problem for the data analysis is that at the higher 
frequencies where lot > 1, the thermally excited quasi- 
particles also contribute to the screening of microwave 
fields and it is then incorrect to use the 1 GHz mea- 
surements of A(T) in order to extract ci(T) from R S (T). 
Effectively, the penetration depth becomes frequency de- 
pendent, a phenomenon that has been observed directly 
in mm-wave measurements on particularly high quality 
thin films |16[j . Ideally this problem can be solved by 
measuring both R S (T) and X S (T) at each measurement 
frequency, but for our higher frequency resonators there 
are considerable technical difficulties with this measure- 
ment. 

It is possible to work around the fact that we only 
have measurements of the real and imaginary part at 
one of the frequencies, as long as there is adequate in- 
formation on the frequency dependence of the real part 
of the surface impedance. The problem is analogous to 
the one commonly faced in far infrared spectroscopy of 
opaque samples, where the reflectance, but not the phase 
of the reflected light, is measured over a wide frequency 
range. Kramers-Kronig relations are the usual solution if 
only one of the optical constants is known, but is known 
over a wide frequency range. Such relations are a conse- 
quence of causality in linear response theory and connect 
one optical constant (such as reflectance) to its partner 
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(phase of reflected light), via an integral over all freque- 
cies. In principle, data at the five microwave frequencies 
presented here could be connected to far infrared mea- 
surements of reflectance in order to perform this analysis 
and extract <J\{uS) and 0*2(0;) . However, such an analysis 
is somewhat difficult because there still exists a substan- 
tial gap in the mm-wave frequency range between our 
highest frequency measurement and the lowest frequency 
far infrared measurements on crystals. Some data has 
been obtained in the mm-wave region using techniques 
such as time domain THz measurements [j^,|l8],|l9| and 
direct infrared absorption |3^] , but these measurements 
have all been performed on films rather than untwinncd 
single crystals. The quasiparticle scattering rate is typ- 
ically much higher in films than it is in single crystals, 
so data taken on such different samples cannot be ana- 
lyzed together in this way. An alternative to a Kramers- 
Kronig analysis is to fit the frequency dependent surface 
resistance to a model, but this is not a very satisfactory 
procedure if one only has data at five microwave frequen- 
cies and one doesn't know the shape of the conductivity 
spectrum a priori. A further problem with both of these 
techniques is that the far infrared data is only available 
at a couple of temperatures below T c , so they can not 
be used to do a complete analysis of the data presented 
here. 

Our main approach to analyzing the surface resistance 
data is to use the R s (T) measurements at 1 GHz to arrive 
at an estimate of how much screening by the thermally 
excited quasiparticles must be included when extracting 
<ti(T) from R S {T) at higher frequencies. Although the 
procedure involves some assumptions about the model 
for the screening, we will show that the corrections are 
small enough that the effect of uncertainty in the choice of 
model does not significantly affect the conductivities that 
we extract in the analysis. We begin by writing down a 
general 2-fluid expression for the microwave conductiv- 
ity that includes contributions to the real and imaginary 
part from both the superfluid and normal fluid (the con- 
ductivity due mainly to thermally excited quasiparticles), 



a(ui,T) = ois 

„2 



n s e 



- «C2S 

.n 3 e 2 



0~1N 
0~1N 



(3.1) 



where cis and a^s are the real and imaginary parts of 
the superfluid conductivity and <j\n and o~2N are the real 
and imaginary parts of the normal fluid conductivity. 
The superfluid contribution consists of a delta-function 
at lu = with an oscillator strangth given by the su- 
perfluid density divided by the effective mass (n s /m*), 
and a Kramer-Kronig related imaginary part that is the 
inductive response of the superfluid at nonzero frequen- 
cies. This superfluid response can be expressed in terms 
of the penetration depth using n s e 2 /m* — (/ioA 2 (T)) _1 , 
and away from to = the delta-function can be neglected, 
leaving 



cr(w,T) = <j 1n (oj,T) - i 



cr 2 Ar(w,r) + 



/z wA 2 (T) 



(3.2) 



Thus, in general the real part of the conductivity comes 
from the normal fluid response and the imaginary part 
has contributions from both the normal and superfluid, 
although the superfluid dominates the imaginary part at 
low frequency. 

In the limit of local electrodynamics the connection 
between this conductivity and the surface impedance is 
made via 



Rs + iX s 



<7l — i<72 



(3.3) 



Thus, at 1 GHz where we have simultaneous measure- 
ments of both R s and X s , o~\ and o~2 can be extracted 
using 



cr 2 = ^qoj 



RsXs 



X 2 



■x 2 ) 2 

R 2 



(3.4) 



(R 2 + X 2 ) 2 



At higher frequencies, where we only have measure- 
ments of Rs(T), it is useful to write a\ in terms of R s 
and (T2 in the following way p9| : 



± 



0-20- 1 



1/2' 



4 



1/2 



(3.5) 



with the sign choice + (— ) for o\ > (<)v / 3o'2 

2RV 



and where a s 



In the normal state, where o\ er 2 at low frequencies, 
this expression reduces to the classical ski n eff ect result 
(Ji = ii u/2R 2 s . The expression in Eq. |3.5| continues 



to be valid right through the superconducting transition 
and can be used to extract <j\ from measurements of R s 
provided that one also has some information on 0-2 ■ The 
extent to which one can use this expression in the super- 
conducting state depends upon the relative importance 
of the normal fluid contribution to cr 2 . In regimes where 
this contribution is small, 02 is simply given by the su- 
perfluid response cj 2 s = (MowA 2 (T)) _1 which can be cal- 
culated from the penetration depth measurements. To 
clearly illustrate the problem that occurs when the nor- 
mal fluid contribution to 02 is not small, we examine a 
particular version of the 2-fluid model where the normal 
fluid conductivity is assumed to have a Drude frequency 
dependence: 



0-\N — 1&2N 



1 + (out) 2 1 + (wr) 2 



(3.6) 



where r is the scattering time of the normal fluid and 
n n /fn* is the normal fluid density over the effective mass. 
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The ratio of the normal fluid and superfluid contributions 
to the effective screening is then 



&2N 
&1S 



M 2 



n s 1 + (uj T f 



(3.7) 



The normal fluid contribution to screening is thus unim- 
portant at low frequency {lot <C 1) or when the normal 
fluid density is small (n n /n s -C 1, which occurs at low 
temperatures) . Conversely, difficulties arise when cot > 1 
and n n /n s is not small, which is likely the case for our 
data at 13 and 22 GHz in the temperature range from 
20-40 K, and for the 75 GHz data from 20-60 K. 

For the conductivities that we will show below we 
take the following approach to estimating the normal 
fluid contribution to 02 . At 1 GHz g\ is calculated di- 
rectly from the simultaneous measurements of R s (T) and 
X S (T). Then, in the spirit of our earlier analysis of this 
type of data, we extract the temperature dependent nor- 
mal fluid scattering time r(T) using the Drude model 
above (Eq. 3.6). The normal fluid density is calculated 
from the penetration depth measurements by assuming 
that 

n„e 2 , , n,e 2 , . n„e 2 , 

-V ( T ) + -V ( ) = "V ( T = °) ■ ( 3 - 8 ) 



which amounts to assuming that all of the low frequency 
oscillator strength is being shifted from the normal fluid 
to the superfluid as the temperature is decreased. The 
t(T) and n n (T) derived from the 1 GHz data can then be 
used to make an estimate of the normal fluid contribution 
to (T2 needed to analyze the data at higher frequencies. 
When this analysis is performed, we find that the normal 
fluid screening influences the extraction of a\ from R s (T) 
only at the level of 20% or less. This can be understood if 
one notices in the R S {T) that lot ~ 1 at temperatures of 
50 K or less at which point the normal fluid density has 
already fallen below 20% of the total oscillator strength. 
Thus, we find that the effect of the normal fluid screening 
is actually fairly small in the microwave range, so any un- 
certainty associated with the choice of a Drude model to 
account for this screening does not have a serious impact 
on the values of <j\ that we derive from the data at 13, 
22, and 75 GHz. Furthermore, the values of r(T) that we 
use in the analysis above turn out to be consistent with 
the conductivity spectra <J\(iS) that will be discussed be- 
low, so the correction for screening by the normal fluid 
is self-consistent. 

Fig. |^ shows the conductivities extracted from the 
R S {T) data of Fig. 0, using the methods described above. 
The sharp upturn at T c marks the presence of super- 
conducting fluctuations, which have been discussed in 
greater detail by Kamal et al. || and Anlage et al. |3|] 
and are not the main focus of the work presented here. A 
number of previous measurements on earlier generations 
oiYBa2Cu^Or-s crystals have shown the presence of an 
extra sample-dependent peak in <Ji(T) just below T c , a 
feature that was discussed by Olson and Koch |3J] and 
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FIG. 4. The temperature dependence of the a-axis mi- 
crowave conductivity of Y BC12CU3O6.993 extracted from the 
surface resistance measurements of Fig. 1. The sharp spike 
near T c is a result of superconducting fluctuations and the 
broad peak at lower temperatures is caused by the increase 
in the scattering time of thermally excited quasiparticles. 



by Glass and Hall |35[ who attributed it to having a sam- 
ple with a broadened transition. A similar feature was 
also reported recently by Srikanth et al. for BaZrO^- 
grown crystals |36|, but since we see no sign of this in 
our new high purity samples, we conclude that such fea- 
tures are associated with a spread in T' c s in the surface of 
the sample. The main feature that we do observe in the 
conductivity is that <J\(T) has a large broad peak, rising 
to nearly 25 times the normal state conductivity. The 
conductivity rises a factor of two higher than was seen in 
measurements at 2 and 4 GHz on an earlier generation of 
crystals grown in yttria-stabilized zirconia crucibles |30"[ . 
The peak also appears at a lower temperature, 25 K in- 
stead of the 35 K turnover observed in the lower purity 
crystals. This effect of very low levels of impurities is 
not consistent with the early suggestion of Klein et al. 
p7| that this feature is the result of BCS-type coherence 
effects. A coherence peak, essentially a density of states 
effect, is observed near T c in <Jx(T) of s-wave supercon- 
ductors such as Pb |||] . However, the peak observed here 
at 1 GHz is much too large and too low in temperature 
to be attributed to such an effect. Furthermore, a strong 
coherence peak has not been seen in NMR measurements 
of T\ in this material |39| nor is it expected in a d-wave 
superconductor, the now widely accepted pairing state of 
YBa 2 Cu 3 7 - S . @0 

Thus, in the absence of strong coherence effects, we 
have attributed the rise in o\ (T) below T c to a rapid in- 
crease in the scattering time r of thermally excited quasi- 
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co(10 s ) 

FIG. 5. The conductivity spectrum at 4 selected tempera- 
tures between 4 and 20 K, extracted from the e>"i(T) curves 
of Fig. 4. In this temperature regime, the conductivity due 
to thermally excited quasiparticles has a nearly temperature 
independent width of 8 GHz and a nearly temperature inde- 
pendent shape that is close to a Drude lineshape (the lines 
are Drude fits). 
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FIG. 6. The conductivity spectrum at 3 selected temper- 
atures above 20 K, extracted from the o"i(T) curves of Fig. 
4. Above 20 K, the width of the conductivity peak broad- 
ens rapidly, stretching out of the microwave frequency range 
above 55 K. These conductivity spectra continue to be rea- 
sonably well fit by Drude lineshapes, as shown by the lines in 
the figure. 



particles, as discussed in the Introduction. The fact that 
this peak rises higher and turns over at a lower tempera- 
ture in the new, higher purity crystals is consistent with 
this interpretation. That is, in the higher purity sample, 
t runs into its impurity limit at a somewhat lower tem- 
perature than it did in the earlier generations of crystals 
and this impurity-limited scattering time is very large in 
the new _BaZr03-grown crystals. An estimate of the in- 
crease can be made as follows. The penetration depth 
measurements indicate that more than 90% of the nor- 
mal fluid density is gone at 25 K (see Fig. |3|), so the 
25-fold increase in o\ between T c and 30 K implies an 
increase in r by at least a factor of 250 over the scatter- 
ing time just above T c . This is actually an underestimate 
because not all of the far infrared oscillator strength in 
the normal state ends up condensed into the superfluid 
at low temperatures. Such a huge increase in lifetime 
is consistent with the relaxation effects observed in the 
data at 13 GHz and higher. As a rough illustration of 
this agreement, far infrared measurements indicate that 
lot ~ 1 at about 3000 GHz just above T c , so a 300-fold 
increase in t below T c would mean lot ~ 1 at about 10 
GHz, leading to the considerable frequency dependence 
in <j\(lo) that we observe in the microwave range. 



IV. CONDUCTIVITY SPECTRA AND 
QUASIPARTICLE LIFETIME 

The evolution of the conductivity with temperature is 
better illustrated in Figs. [| and p| where we show the 
conductivity spectrum <j\ (lo) at several different temper- 
atures. In fact, the central technical achievement of this 
work is that we now have measurements at enough fre- 
quencies that both the shape of <Ji{ui) and its tempera- 
ture dependence are quite clear. Fig. |5| shows the con- 
ductivity spectrum at four different temperatures below 
25 K, the temperature range where we have argued above 
that the conductivity is dominated by thermally excited 
quasiparticles that are being scattered by a low level of 
impurities. We will not concentrate here on the lowest 
temperature data, where the loss becomes very small and 
difficult to measure accurately by most cavity perturba- 
tion techniques. Fig. || shows that from 4 to 20 K, the 
conductivity consists of a very narrow peak whose width 
is about 8 GHz. The lines in the figure are fits to a 
Drude lineshape (the real part of Eq. 3.6), demonstrat- 
ing that the conductivity spectrum of the thermally ex- 
cited quasiparticles is Drude-like. There are deviations 
from this lineshape; in particular, the data at the low- 
est frequencies tends to lie above any Drude curve that 
manages to fit the data at 13 GHz and above. That is, 
cri(u>) seems slightly more cusp-like than a Lorentzian 
curve. Despite these somewhat subtle details, the nearly 
Lorentzian lineshape largely confirms the ansatz origi- 
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FIG. 7. A comparison of the normal fluid oscillator 
strength determined in two ways; from Drude fits to spec- 
tra like those of Figs. 5 and 6 and from the disappearance of 
oscillator strength in the superfluid response, as measured by 
1/A 2 (Fig. 3). 



nally used by Bonn et al. |29| to analyze the early R S (T) 
measurements. Now, however, it is possible to glean an 
effective scattering rate 1/r by fitting <Ji(oj) to a Drude 
lineshape, rather than assuming one. The fits shown in 
Fig. |^ all yield a lincwidth close to 8 GHz, so the dom- 
inant change in the spectra in this temperature range 
seems to be an increase in oscillator strength, due to the 
shift of spectral weight from the superfluid response (a 
^-function at lu=0) to the microwave conductivity. 

One sees in Fig. ^ that above 25 K the conductivity 
peak starts to broaden rapidly and by 60 K the width 
becomes much greater than the frequency range of the 
microwave measurements, so we have no direct measure 
of the width and shape of <Ji(u) above 60 K. Again, in 
this figure the lines are fits to a Drude lineshape. The 
fits are not perfect, but the flattening of the curve at the 
lowest two frequencies and the rapid roll-off with increas- 
ing frequency are key features of a Drude conductivity 
spectrum. Although the Drude lineshape is not a perfect 
fit to the data (even with the 5-10% uncertainty in the 
absolute value at the various frequencies), overall the fits 
between 4 and 55 K are good enough to consider plotting 
the fit parameters as a function of temperature. 

An interesting cross-check of these fits to a\(uo) is 
to compare the oscillator strength in the normal fluid 
conductivity peak, which is given by the fit parameter 
n„e 2 /m*, to the superfluid density n s e 2 /m* = (^ D A 2 ) -1 
extracted from the penetration depth measurements. If 
one assumes that all of the oscillator strength ends up 
in the superfluid 5-funtion as T— >0, then the superfluid 
density is related to the normal fluid density via Equa- 




T(K) 

FIG. 8. The scattering rate of the thermally excited quasi- 
particles, as inferred from the width of Drude fits to the con- 
ductivity spectra of Figs. 5 and 6. 



3.£ . Fig. shows a comparison between the normal 



tion 

fluid density inferred from the fits to the ai(w) peaks 
and the normal fluid density inferred from the penetra- 
tion depth via Eq. 3.8 . This figure indicates that the 



normal fluid density does track the decline of the super- 
fluid density with increasing temperature, which lends 
some degree of confidence to this parameterization of the 
data. The agreement at low temperatures is excellent, 
with an increasingly serious deviation betweeen the two 
above 30 K. The lack of agreement at the higher tem- 
peratures is an indication that the Drude lineshapes do 
not completely keep track of where all of the oscillator 
strength is going as temperature increases. There is some 
deviation at low temperatures, taking the form of a nor- 
mal fluid oscillator strength that is extrapolating linearly 
towards a non-zero value as T— >0. This is an indication 
of the presence of residual conductivity in the low tem- 
perature, low frequency limit, which is expected for a 
d-wave superconductor p2| . 

Although the comparison of the oscillator strengths 
shown in Fig. [?] indicates that the Drude fits to ai(u>) do 
not perfectly account for the redistribution of oscillator 
strength, the shape does provide a reasonable measure of 
the width of the peaks from 4 to 55 K. The temperature 
dependent width coming from these fits, which we have 
previously interpreted as the scattering rate of the ther- 
mally excited quasiparticles, is shown in Fig. |[ One of 
the key results of these measurements is that the width 
of the normal fluid peak is very small, only about 8 GHz, 
and it is nearly temperature-independent up to 20 K. The 
8 GHz width suggests a very long quasiparticle scattering 
time of 2xl0 _11 s. In previous work performed on lower 
purity samples measured at only 2 frequencies (4 and 
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35 GHz) a similar conclusion was drawn regarding the 
spectrum of the thermally excited quasiparticles. In that 
data it was suggested that the width was about 30 GHz, 
with little temperature dependence of the scattering rate 
up to somewhat higher temperature, about 30 K. The 
narrower width of 8 GHz lasting only up to about 20 K 
is consistent with the higher purity of the new BaZrO^- 
grown crystals. That is, the lower impurity limit to the 
width makes the broadening due to inelastic scattering 
apparent at somewhat lower temperatures. 

Above 20 K the scattering rate increases very rapidly, 
broadening the peak in a±(uj) out of the microwave range 
by 60 K. Eventually, this width will broaden out into the 
far infrared conductivity spectrum in the normal state, 
a width of ~ 2ksT which is about 3600 GHz just above 
T c . That is, the quasiparticle scattering rate falls by a 
factor of about 500 from T c down to the impurity limit 
encountered at 20 K. Thus, the exceptionally low level 
of defects in the new BaZrO^ grown crystals, coupled 
with the measurements at 5 microwave frequencies, al- 
lows us to determine the temperature dependence of the 
inelastic scattering rate over a wide temperature range. 
Qualitatively similar results have been obtained in mea- 
surements on thin films by high frequency microwave and 
THz techniques ||[l9]]. However, the much higher level of 
defects in such samples limits the temperature range over 
which the evolution of the inelatic scattering rate can be 
observed. In the high purity crystals, the temperature 
dependent scattering can be tracked all the way down to 
20 K. 



V. DISCUSSION 

For the purposes of discussing the data, we have di- 
vided the conductivity spectra up into the two regimes 
discussed in the previous section. In the range below 
about 20 K where the width of the peak in a\(uo) is 
narrow and nearly temperature independent, studies of 
samples over a wide range of purities indicate that this 
regime is governed by thermally excited quasiparticles 
being scattered by impurities or other defects p0| , p0[ . 
Waldram et al. have pointed out the possibility that 
non-local effects might come into play in the conductiv- 
ity in this regime, leading to an effective scattering rate 
that is not influenced by the density of residual impurities 
p3| . However, the considerable narrowing of <Jx(u>) that 
we have observed upon going from YSZ-grown crystals 
to the higher purity BaZrO^-growii crystals indicates 
that the samples are still in a regime where impurities 
play a role in the low frequency scattering. We have 
previously suggested an intuitively appealing way to in- 
terpret the spectra in this regime, based on a specific 
version of the two fluid model p9[ . In this phenomeno- 
logical picture, the quasiparticles excited near the nodes 
in the energy gap have a temperature-independent scat- 
tering rate due to elastic scattering by impurities and 



a conductivity spectrum with a Drude lineshape whose 
width is set by this scattering rate just like impu- 
rity scattering in a normal metal. In the high purity 
samples this would correspond to a strikingly long 
mean free path of 4 \xm if one takes the Fermi velocity 
to be up = 2 x lQ 7 cm/s. In this particular two fluid 
model, the only source of temperature dependence in the 
low temperature microwave conductivity is the density 
of the thermally excited quasiparticles, which increases 
linearly with temperature. This is a straightforward con- 
sequence of the linear dispersion of the gap function near 
the nodes and is also intimately connected to the linear 
temperature dependence of the penetration depth (they 
are Kramers-Kronig related) . 

Parameterizing the normal fluid response in this way, 
with a temperature dependent density and a scattering 
rate, has been partly justified by calculations of the mi- 
crowave conductivity of a d x 2_ y 2 superconductor [ Ell . 
However, it is not so obvious that a temperature indepen- 
dent scattering rate is expected for impurity scattering 
of these thermally excited quasiparticles near the nodes, 
since they obviously differ greatly from free carriers at 
an ungapped Fermi surface. In particular, it has been 
pointed out by Hirschfeld et al. that elastic impurity 
scattering in this situation should lead to a frequency 
and temperature dependent scattering rate because of 
the restricted phase space into which the quasiparticles 
at the nodes can scatter p4] , ^6| . So, with this possible 
conflict between theory and the phenomenological model 
in mind, we will make some more detailed comparison be- 
tween the microwave conductivity data and the relevant 
theoretical calculations. 

The transport properties (both microwave conductiv- 
ity and thermal conductivity) of a d x i_ y i superconductor 
have been the subject of considerable theoretical effort 
recently [EJ|,|4| Q . This work builds on earlier calcu- 
lations of the transport properties of anisotropic super- 
conductors, aimed primarily at explaining and predicting 
the properties of heavy fermion superconductors |^l| , ^2| . 
The high temperature superconductors now offer an op- 
portunity to test these calculations in a situation where 
we seem to have identified a relatively simple anisotropic 
pairing state. Such comparisons are somewhat complex 
because the question of transport properties at low tem- 
peratures is inherently a question of understanding im- 
purity effects. This is especially the case for anisotropic 
superconductors because the presence of impurities has a 
strong impact on the excitation spectrum near gap nodes, 
particularly in the limit of unitary scattering. 

A key effect of impurities in an anisotropic supercon- 
ductor is to produce a band of impurity states with some 
width 7, thus giving the superconductor a non-zero den- 
sity of states at the Fermi level. One surprising conse- 
quence of these states is a universal conductivity limit 
at low frequency as T — > 0, first pointed out by PA. 
Lee [EJ|. This residual conductivity is independent of 
impurity concentration, the result of a cancellation be- 
tween the density of states induced by the presence of 







the impurities and the transport lifetime associated with 
those states. A version of this universal limit has been 
observed by Taillefer et al. in thermal conductivity mea- 
surements of pure and Zn-doped YBa^Cu^O^.^ below 1 
Kelvin [^3|. To the best of our knowledge, this limit has E 
not yet been definitively observed in microwave conduc- 'qj 
tivity measurements, in part due to sensitivity problems r o 
in the type of cavity perturbation measurement being ^ 
discussed in this article. Instead, our main concern here ® 
will be the behaviour of this conductivity as temperature 
and frequency are increased, which involves conductivi- 
ties that are substantially larger than the T— >0 limit and 
are thus easily measurable with the methods discussed 
above. This microwave conductivity has been the subject 
of considerable theoretical effort, including both numer- 
ical work and a number of analytical results in certain 
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One particularly well studied limit involves the electro- 
dynamic properties at temperatures and frequencies be- 
low the impurity bandwidth 7, in the unitary scattering 
limit (scattering phase shift — > 7r/2). In this limit the im- 
purity bandwidth is given roughly by 7 ~ (TA) 1 / 2 where 
A is the magnitude of the gap and T is the elastic scat- 
tering rate that the impurities would contribute to the 
normal state resistivity. For w,T < 7, where the trans- 
port properties are dominated by this impurity band, it 
has been shown that both o\ and A vary as T 2 Q,^5|. 
This quadratic behaviour has been seen in Zn-doped sam- 
ples of YBC12CU3O6.95, where it was found that at a Zn 
impurity concentration as low as 0.15 % the crossover 
energy scale is already 7 > 4 K jl9|,[30|]. However, Zn 
substitution for planar Cu's is the only impurity that we 
have found that clearly gives this unitary scattering be- 
haviour. Ni substitution for Cu, Ca substitution for Y, 
and the chain oxygen vacancies all seem to have much 
weaker effects, even at defect levels of 1% or more 
The previous generation of YSZ-grown crystals showed 
only slight curvature in A(T) and o~\{T) below 4 Kelvin 
and the new BaZr03-grown crystals show no hint of T 2 
temperature dependence down to 1.2 Kelvin. The rela- 
tive rarity of this quadratic behaviour (though it is com- 
mon in thin films for reasons that remain unclear 54 1) 
leads us to consider the opposite limit for the strength of 
the scattering, the Born limit. 

For impurity scattering in the Born limit, the crossover 
energy scale is exponentially small, so one does not nec- 
essarily expect to see the universal conductivity limit un- 
til measurements are performed well below 1 Kelvin. In 
fact, in the microwave measurements presented here we 
are not necassarily at low enough temperature or fre- 
quency to observe any simple limiting behaviour. So, 
in this situation we compare our results qualitatively to 
numerical calculations performed in the Born limit by 
Hirschfeld et al. |iq| . They found that at very low fre- 
quency o~\(T) rises rapidly from the universal zero tem- 
perature limit to a much larger conductivity that depends 
upon the impurity scattering rate. It then remains fairly 
temperature independent until inelastic scattering pro- 
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FIG. 9. This detailed view of the temperature dependence 
of the microwave conductivity below 20 K shows a gradual 
evolution of the shape of of ai(T), from concave up at high fre- 
quencies to concave down at the lowest frequencies. The quite 
linear temperature dependence seen here at 13 GHz seems to 
be an intermediate behaviour. 



cesses become important. At higher frequencies 01 (T) 
becomes smaller and moves through a whole range of be- 
haviours, varying from mostly sub-linear in T at low fre- 
quency, through a quasi-linear temperature dependence 
at intermediate frequencies, to a faster than linear tem- 
perature dependence at high frequencies. Figure ^ shows 
behaviour in the measured microwave conductivity that 
is similar to this Born limit result in some of its qualita- 
tive features. The overall conductivity has a magnitude 
that varies with purity, from the quite high values seen 
here, through to very low, flat conductivities observed 
in Ni-doped samples ]3(J. More importantly, here we do 
clearly see for the first time that the linear behaviour 
of <7i (T) is an intermediate behaviour, albeit one that 
survives over a substantial range of frequency and tem- 
perature. The evolution in the shape of o~i(T) at low 
T clearly falls outside of our phenomenological model, 
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which would have predicted a linear temperature depen- 
dence at all of the frequencies shown in Fig. ||. That 
is, if <7i (u>) were really Drude-shaped with a temperature 
independent width, then <J\{T) would exhibit the same 
temperature dependence at all frequencies; namely, the 
linear temperature dependence of the normal fluid den- 
sity. 

At our lowest frequencies the data does move towards 
sub-linear as expected in the Born limit. However, one 
perhaps important difference from the theoretical calcu- 
lations is that the trend does not seem to continue below 
2 GHz. The data stops evolving towards the expected 
low frequency behaviour, which is a rapid leap upwards 
to a constant value. The reason for this is not yet clear, 
so we are left for the moment with certain qualitative 
features that seem in accord with Born limit scattering. 
Other features of the data seem to echo the expectations 
of a conductivity spectrum with a Drude-like lineshape 
that has a temperature independent width; namely, the 
similarity in the 1 and 2 GHz curves and the relatively 
large frequency and temperature range over which <j\ (T) 
seems to be linear in T. It is the latter aspects of the data 
that lead to the phenomenological model being a fairly 
good description of the conductivity spectra. 

Although the foregoing discussion indicates that the 
<7i(tL>) spectra are not quite Drude shaped, there is still a 
characteristic width to the peaks that is reasonably well 
parameterized by the Drude fits. Thus, the plot of scat- 
tering rates shown in Fig. || still provides a reasonable 
measure of the narrowness of the peak at low temper- 
ature and its rapid broadening above 20 K. In Fig. [h] 
we plot this width, which we identify with the quasipar- 
ticle scattering rate, versus T A . Just above 20 K, the 
initial onset of this increase in scattering appears to be 
at least as rapid as T 4 and rises even more quickly at 
higher temperatures. This type of rapid temperature de- 
pendence of the quasiparticle scattering time would be 
expected in any situation where the inelatic scattering 
comes from interactions that become gapped below T c . 
A number of early calculations tackled the problem of the 
collapse of the scattering rate below T c in this way. Early 
on, Nuss et al. explained the peak in their THz conduc- 
tivity measurements in this manner Littlewood and 
Varma studied this type of effect in a marginal Fermi 
liquid |5jJ, the idea being that below T c a gap opens 
up in the scattering spectrum. Statt and Griffin simi- 
larly studied the effect of the opening of a gap in the 
spectrum of spin fluctuations |5q| . All of this work pre- 
dated the solid identification of the d x 2_ y 2 pairing state, 
so isotropic s-wave gaps were assumed in the calculations. 
Quinlan ct al. studied a model of quasiparticle scattering 
in which the lifetimes were associated with spin fluctua- 
tion scattering. They studied the effects of both s— wave 
and d— wave gaps opening up in the spin fluctuation spec- 
trum |57]] . Since we now know that the gap in this mate- 
rial has d x 2_ y 2 symmetry, this latter calculation is most 
directly relevant to our measurements here. In particu- 
lar they found that at temperatures well below T c , the 




T (10 K ) 

FIG. 10. When the temperature dependence of the scatter- 
ing rate of thermally excited quasiparticles is plotted against 
T 4 as shown above it is clear that the onset of inelastic scat- 
tering is very rapid, at least T 4 and even faster above 30 K. 

quasiparticle lifetime increases as T 3 and even faster than 
this as T c is approached. In a comparable temperature 
range, the scattering rate that we extract from the width 
of <t\{uj) is closer to T 4 . Thus, the temperature depen- 
dence of the inelastic scattering rate seems to be about 
one power of T faster than the lifetime calculations based 
on a gapping of the spin fluctuation spectrum. 



VI. CONCLUSIONS 

We have presented here the most complete microwave 
measurements yet obtained on a crystal of a high tem- 
perature superconductor. These spectra reveal a wealth 
of detail regarding the conductivity spectrum ai(u>) of 
the thermally excited quaiparticles in the superconduct- 
ing state. We find that ai(u>) has a Drude-like shape. 
Although in detail there are deviations from this shape, 
there is nevertheless a well defined characteristic width 
which can be associated with the scattering time of the 
quasiparticles. Below T c the width collapses rapidly and 
it becomes easily discernible in the microwave spectral 
range at temperatures below 55 K. In the range between 
55 K and 20 K we find that the collapse of the scattering 
rate varies at least as fast as T 4 . By 20 K the width 
becomes extremely narrow and nearly temperature inde- 
pendent. This narrow width of only 8 GHz corresponds 
to a mean free path as high as 4 /im if we interpret the 
width as being a direct measure of the elastic scatter- 
ing rate due to impurities. We find that some features of 
crj (u>, T) below 20 K are in accord with quasiparticle scat- 
tering in the Born limit for a d x 2_ y 2 superconductor, in 
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particular, a gradual evolution of the shape of <J\{T) from 
sublincar T dependence at low lu, to quasilincar and then 
faster at higher u>. However, there remain discrepancies 
that might best be settled by a detailed numerical calcu- 
lation aimed at fitting the observations presented here. 
Such fits must come to grips with the observation that 
the behaviour of o\ (uj, T) below 20 K is reasonably well 
described by a model involving quasiparticle scattering 
with a temperature-independent lifetime. 
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